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NOMENCLATURE 

analytical limiter input (used in describing function analysis) 

analytical limiter value (used in describing function analysis) 

describing function gain 

deadband of rudder pedal deflection 

deadband in rudder pedal integrator 

rudder pedal to rudder gearing function and limiter 

digital computer sampling frequency 

full authority series servo gain 

open loop transfer function 

moment of inertia about aircraft z-axis 

root locus gain 

control mode phase-out gain 

gain of roll attitude feedback into yaw controller 
parallel/series servo mode combined analysis gain 

■ - So “ Sc ^ ^ 

parallel servo gain 

parallel/series servo mode combined analysis gain * (K x K..) 

SS 

series servo gain 
actuation dynamics torque gain 
additional parallel servo gain 
yaw attitude feedback gain 
yaw rate feedback gain 
pedal Integrator gain 
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^ combined analysis gain = (K^ x 

combined analysis gai, = (Kiq + ^20) 

K3 full authority series servo mode combined analysis gain 

= (K30 X * 1 X y FCN X K.^oRQUe/^e) 

parallel/series servo mode combined analysis gain 

“ ^*^TORQUe/^z^ 

K]o controller forward gain component 

K20 controller forward gain component 

Y 0 y^w controller coupling gain 

LMry position limit of rudder pedals 

yaw controller coupling limit 
SRFIMF state-rate-feedback- implicit-model-following 

T digital computer sampling period 

3 aircraft sideslip angle 

time derivative of inertial sideslip angle 
6 x^ pilot control input to yaw controller (pedal Input) 

actuation model time constant 
flight controller compensation time constant 
aircraft roll attitude 
i(; aircraft yaw attitude (heading) 

^ aircraft yaw rate 

aircraft yaw acceleration 
experimental limiter output 
experimental limiter input 
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AN ANALYSIS OF A NONLINEAR INSTABILITY IN THE IMPLEMENTATION 


OF A VTOL CONTROL SYSTEM DURING HOVER 
Jeanine M. Weber 
NASA-Ames Research Center 


SUMMARY 


An analysis has been conducted to determine the contributions to non- 
linear behavior and unstable response of the model following yaw control 
system of a VTOL aircraft during hover. The system was designed as a state 
rate feedback implicit model follower that provided yaw rate command/heading 
hold capability and used combined full authority parallel and limited 
authority series servo actuators to generate an input to the yaw reaction 
control system of the aircraft. Involved in the analysis were linear and non- 
linear system models, and describing function linearization techniques to 
determine the influence of input magnitude and bandwidth, series se* vo 
authority and system bandwidth on the control system instability. Results 
of the analysis describe stability boundaries as a function of these system 
design characteristics. 


INTRODUCTION 


The use of advanced control systems that provide stabilization and 
command augmentation for attitude, translational velocity, and position 
control have been shown to reduce pilot workload for VTOL hover operations 
(references 1 and 2). One such example is the state-rate-feedback-implicit- 
model-following (SRFIMF) concept examined in reference 3, which achieves 
model-following fidelity through the feedback of acceleration, rate and atti- 
tude signals* In the ground simulation experiment of reference 3, this 
concept was shown to enhance the capability of the pilot-aircraft system to 
perform a demanding hover (and deceleration) task. 

In-flight evaluation of several such control systems, including SRFIMF, 
with a VTOL research aircraft has been proposed* Toward that end, it is 
anticipated that the control system of an existing operational VTOL aircraft 
will be modified to permit incorporation of these concepts* In order to 
minimize the cost of implementation, a simplex electronic control system with 
manual override capability has been suggested. Safety considerations would 
be satisfied with a parallel/series servo arrangement of the actuators such 
that ea^h servo is limited to contain a runaway failure. The parallel servo 
would be full authority and have a limited rate-of-actuation* The series 
servo would be position limited and capable of a high rate-of-actuatlon. 



The reference 3 ground simulation assumed a full-authority fly-by-wire 
control system and did not address possible influences of the parallel/series 
servo mechanization on the SRFIMF control system. Hover simulation of the 
SRFIMF controller implemented in this manner during a control system failure 
study at Ames Research Center in August, 1979 (reference 4), revealed a non- 
linear instability in the yaw axis. An electrical limiter in the control 
system, designed to limit the servo input upstream, s.iturated and cycled 
causing the divergent oscillation of the output signals. Although during the 
simulation, the instability was apparent only in the yaw axis, it is antici- 
pated that there is a potential problem in any channel using the parallel/ 
series mechanization. 

The purpose of this study was to examine the causes and characteristics 
of this nonlinear instability. A describing function analysis technique 
(reference 5) was used in the study. Simplifies'’’ linear and non-linear models 
of the controller, actuation and airframe dynamics were used to isolate and 
analyze the problem experimentally. The analysis encompassed only the wings- 
level hover flight condition. 

This report presents the results of this analysis, including a descrip- 
tion of the model with the appropriate simplifications, and analyses of the 
linear and nonlinear systems. The Nonlinear System Analysis portion Involves 
a description of the system, theoretical and experimental approaches to under- 
standing nominal system response, and experimentally establishing stability 
boundaries by varying system configuration. 


MODEL DESCRIPTION 


The SRFIMF yaw control system has been simplified for purposes of 
analysis. Included in the simplified model are the yaw flight controller, 
the option for selecting either parallel/series or full authority series 
servo implementation, actuation dynamics and yaw response of the aircraft 
which neglects aerodynamic effects. 

Figure 1 shows the SRFIMF yaw control system as taken from figure 14 of 
reference 3. Shown in figure 2 is the yaw control system used in the 
analysis, which includes the following simplifications from the figure 1 
model: deletion of the electrical deadbands on the pedal Inputs to the sys- 
tem (DB(i;u dnd deletion of the sideslip command mode (located down- 
stream or the Bx ^ inputs) which is active only for flight at or 

above 30 knots, and deletion of the tan^ input path which is inapplicable to 
this wings-level analysis* 

Figure 3 shows a block diagram of the simplified overall control system 
including servo Implementation and actuation dynamics. The implementation 
flag selects either the parallel/series servo mode (1) or the full authority 
series servo mode (0). The limiter shown in figure 2 limits the servo 

input signal in the case of a runaway failure. A digital representation of 
the model shown In figures 2 and 3 was used in the analysis. 
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The problem has been analyzed in two parts; Initially tlie limiters were 
removed and the linear response was studied in each of the parallel/series 
and full authority series servo modes. The p >se of this was to understand 
the linear system response, and to verify the .aplified Fortran model. The 
limiters were then introduced and the nonlinf response was analyzed theo- 
retically with the describing function technique and experimentally using 
the time histories generated by the Fortran program. 


LINEAR SYSTEM ANALYSTS 


Linear System Classical Stability Analysis 

This section describes the response characteristics of the linear system 
in terms of classical analysis techniques. 

With further simplification, the block diagrams in figures 2 and 3 have 
been reduced to figures 4 and 5, representing the parallel/series and full 
authority series modes respectively. 

Gain and signal equivalences between SRFIMF and the analysis, as well as 
the nominal gain values are given in the Appendix. 

From figures 4 and 5 the open loop transfer functions of each mode can 
be written: 


GH(s) 


Parallel/Serles 
Servo Mode 
(Linear System) 





GH(s) 


Series Servo 
Mode 

(Linear System) 



( 1 ) 


( 2 ) 


where K Is the open loop gain factor. 

A root locus of the parallel/series mode is shown in figure 6 with 
partial root loci of each mode in figures 7 and 8 showing only the Imaginary 
axis region. From figures 7 and 8 it is apparent that the system is condi- 
tionally stable in each mode, thus gain reduction is sufficient to destabi- 
lize the response. Note that in the parallel/series mode, the high frequency 
eigenvalues cross the imaginary axis at a higher gain value than in the series 
servo mode, thus the series servo mode remains stable at a lower effective 
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system gain. From equations (1) and (2) it can be seen that the series servo 
mode has three open loop integrators, whereas the parallel/series mode has 
four. In the nonlinear system, the extra integrator associated with the 
parallel/series servo implementation introduced a further destabilizing 
effect to the system. 

Relative degrees of conditional stability between the modes are also 
illustrated using frequency response techniques. Figures 9 and 10 are Bode 
diagrams of the linear system. The higher gain and phase margin of the 
series servo mode indicate the increased resilience of that mode to system 
gain and phase angle reduction. From the figures it can be seen that the 
series servo mode becomes neutrally stable at the phase crossover frequency 
(2.2 rad/sec); above this frequency the system has a stable response. Due to 
the conditional stabiltiy of the system, if the Bode gain is reduced 21 db 
below the total loop gain of the design, the system becomes unstable. By 
contrast, the parallel/series servo mode has a phase crossover frequency of 
4,5 rad/sec and the total loop gain can only be reduced by 11,6 db without 
destabilizing the system. 


Computer Model Verification 

The nominal sampling period, T, of the full Harrier simulation on the 
Xerox Sigma 9 computer is .05 sec, corresponding to a sampling frequency, fs, 
of 20 Hz. However, when the linearized model (limiters removed) is run at 
T - .05 sec, the characteristics of the time response are not equivalent to 
those predicted by the linear root locus analysis. When the sampling fre- 
quency is Increased by a f^<ctor of ten (T * .005, f^ * 200 Hz), the time 
history characteristics cotrp^jre well to those predicted by the linear analy- 
sis in the region of marginal to neutral stability. Figure 11 shows a com- 
parison of the marginally stable mode characteristics as extracted from the 
time history responses run at T * .05, T * .005 and the linear analysis. 

It is apparent from this comparison that the slower sampling frequency 
reduces the damping of this mode substantially and can lead to erroneous 
conclusions regarding system stability. Thus, for accuracy purposes, the 
analysis of this problem has been accomplished with a sampling period of 
T » .005 sec. 


NONLINEAR SYSTEM ANALYSIS 


Nonlinear Control System Elements 

The non-linear control system Includes limiters on series servo position 
and on the yaw reaction control. 

The purpose of the first limiter is to arrest the aeries servo input to 
the yaw control actuator in the case of a series runaway failure. The limits 
are determined by the percentage of total control power that is acceptab.le to 
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be allotted to the series servo. The practical limits on series servo control 
power are a compromise betveen that required for effective control augmenta- 
tion and that acceptable based on hard-over servo failure. With inputs of 
sufficient magnitudes and above specific frequencies the series servo was 
found to saturate and limit cycle, causing a divergent oscillation of the yaw 
response of the aircraft. 

The yaw reaction control limit simply represents the total yawing moment 
control authority. This authority limit may be reached if the parallel 
actuator is driven to sufficiently large magnitudes by the pilot's control 
Inputs or by unstable response in the yaw axis. 


Describing Function Analysis 

One tool for linearization and analysis of nonlinear control system 
elements is the describing function which is discussed in reference 5. The 
describing function of a simple limiter is: 


where 




A 



a 


gain which replaces limiter in the linearization 
magnitude of signal input to limiter 
value to which output is limited 


(3) 


It can be seen from equation (3) that as the ratio of limiter size to signal 
magnitude, a/A. varies between 1 and 0, the describing function gain, bi/A, 
varies between 1 and 0 also. 


Shown in figure 12 is the parallel/series mode block diagram including 
the nonlinear element and, in equation (A), the open-loop transfer function 
with the equivalent linearized element. (A similar transfer function may be 
written for the series servo mode.) 


GH(s) 


Parallel/Series 
Servo Mode 
(Linearized 
System) 



(A) 
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As Llie limit on servo command is reached and b| /A decre*ises from the 
linear value of I, the effects on the open loop transfer function are: 

1. a reduction in system gain 

2. movement of a pole from the origin along the negative real is. 

By solving the characteristic equation at discrete bj /A values^ the cor- 
responding variation of the closed-loop roots of the system has been deter- 
mined. Figures 13 and 14 show these root locations in each of the paralKl/ 
series and series servo modes respectively. In these figures, it may be seen 
that as bj/A decreases from 1, a low frequency mode is introduced. The 
damping of this mode is reduced progressively as bj/A approaches zero. In 
the parallel/series mode (figure 13), as bj /A approaches zero, the closed- 
loop roots of the system approach those of the open-loop case that has three 
zeros at the origin. This leads to the low frequency eigenvalues being 
forced into the right half-plane. In figure 14, however, the closed-loop 
roots of the series system approach only two zeros at the origin as bj /A 
approaches zero; thus, the root locus branches approach vertically and the 
eigenvalues always remain in the left half-plane. 

It may be concluded that describing function gain value has a signifi- 
cant effect on system response and stability. Given system input, the 

describing function gain value, bj /A, may be found with the transfer function 
between system and limiter input, A/6j^, given in equation (5). 



The complexity of the equation as a function of A prohibited solving 
explicitly for A in terms of Instead, an iterative procedure was 

used which consisted of: 

1. guessing A 

2. calculating the corresponding /A 

3. calculating A « x RHS 

4. comparing calculated A with guess and. If they %rere not equal, 
returning to 1. 
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This iterative procedure was followed untix the guessed value of A from 
step //I converged upcn the calculated value of A from step //3, when the 
corresponding describing function gain value, bj /A, for the given system 
input, was able to be found from equation (3). Note equations (3) and 

( 5 ) show chat although the describing function gain is inherently only magni- 
tude dependent, its use in the closed loop analysis causes the limiter input 
magnitude. A, and hence b^ /A to be irequency dependent also. Thus, the fre- 
quency as well as the amplitude component of system input, < 5 x 4 ^* affected the 
iteration on bj /A. 

It is possible from this relationship to determine Che effect on system 
stability of th^ amplitude and frequency of the p..lots input ^ 6 x 4 ^)* the servo 
command limit (a), and the control system bandwidth as determined by the 
overall system loop gain. 


Effect of Implementation on Time Response 

In this section, the differences in nonlinear dynamic response between 
the parallel/series and series servo implementations are illustrated with 
sample time responses and comparisons of those responses with the theoretical 
predictions are made as a means of assessing the validity of the theoretical 
model. 

From the discussion in the previous section it may be anticipated that 
the discernible effect of parallel/series servo Implementation is a divergent 
response once the magnitude of the pilot's input has reached a sufficient 
level. It may also be anticipated that, although the damping of the series 
servo mode response is greatly reduced, the response remains stable for all 
inputs. 

Shown in figures 15 through 22 is an ex* erlmental comparison between the 
responses of the parallel/series and series only modes at various input 
magnitudes and frequencies for the nominal limiter authority of ±.12 inches. 
Among the time history traces shown are (top to bottom) limiter input, 
limiter output » yaw acceleration, rate and heading and the pilot's pedal 
input. Table 1 summaiizes the theoretically predicted response (l.e.» limilctr 
input magnitude and nonlinear mode frequency and damping) and the results of 
the experimental comparison at each system configuration. 

From figures 15 and 16, it may be seen that the parallel/series system 
is stable with commanded output, ji, in response to each of the low frequency 
inputs. In figure 15, the low amplitude input of .1 in. does not saturate 
the limiter. With an input amplitude of .4 in., figure 16 shows the initial 
response of United to .12 in. with respect to however, the 

response after this transient remains unsaturated. 

Each of these cases was predicted very veil by theory as shown in 
Table 1, the unsaturated responses were predicted by a describing function 
gain of 1.0 (limiter unsaturated), and no destabilising nonlinear mode. 
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The steady state limiter input magnitudes, A or also corroborated 

very well. 

Increasing the frec^uency of the low amplitude input to 8 rad/sec in 
figure 17 shows steady state limiter input, to have increased to .23 in. 

(compared with a predicted value of .24 in.) and i/z/vL ^ stable limit 
cycle. The describing function gain of .56 predicted the introduction of the 
lightly damped, low frequency mode overlaying the system output, i//. From the 
table, the gap between theory and experiment has widened here, with the non-’ 
linear mode predicted to be of higher frequency and higher damping than the 
experimental results shown. In ' -gure 18, increasing the input amplitude to 
.4 in. yielded an unstable limit cycle, 4 >aL» ^ divergent yaw attitude 

response, (r, » -.081). Theory again predicted the initial limiter input 
accurately, but predicted a higher frequency and more highly damped nonlinear 
mode (c * 0). 

To summarize the parallel/series configuration response, it was found 
that an unstable system response may be obtained by increasing either the 
input amplitude or frequency above a sufficient level. Predictions made 
usln/, the describing function linearization technique were very accurate when 
calculating lirriter input magnitude; however, onl 'near mode was con- 

sistently less stable than that predicted by theory. 

The series servo configuration tolerates a higher limiter input signal 
before saturation since the limits are proportional to series servo authority 
(full authority for this configuration). System responses to each of the low 
frequency inputs, shown in figures 19 and 20, snow no saturation and no non- 
linear mode overlaying ;he output. From Table 1 it may be seen that these 
responses were well predicted by theory. Increasing input frequency to 
8 rad/sec in the low amplitude case, figure 21, produces no saturation due tc 
the larger limiter tolerance. (Recall the stable limit cycle in figure 17 for 
a .1 sin (8t) input into the parallel/series mode.) Again, no destabilizing 
nonlinear mode is present in the tine history, consistent with the theoretical 
prediction. Increasing the input amplitude to .4 in. (figure 22) produces a 
stable limit cycle, and a fairly well damped nonlinear mode ■ .54). 

Again, theory predicted a more stable nonlinear mode (C ■ *75) tfuin was found 
from experiment. Recall that the response of the parallel/series system to 
dll input of .4 sin (8t) (figure 18) was an unstable limit cycle, and divergent 
aircraft yaw response. 

Theoretical predictions of series servo mode response were found to be 
approximately as accurate as those of the parallel/series mode. In each case 
limiter input magnitude was predicted very accurately; however, the nonlinear 
mode damping was overpredicted, indicating a less stable system experimentally 
than would be anticipated from theory. 

In this section, it has been found that the describing function linear- 
ization provides a reasonably accurate pre^iiction uf aystem behavior* By 
substantiating the results of the previous tn^ioxetlcal dlscassion, it may be 
concluded that the full authority series configuration, while aff^ted by the 
low frequency nonlinear mode, will always remain stable. Prediction of the 

8 


i 




stability of the parallel/series mode provides the motivation for the balance 
of this report. Recall that the critical differences in system response 
between the implementations were due to: 

1. a lower limiter input magnitude tolerance in the parallel/series 
system, and 

2. the extra open loop pole at the origin introduced by the 
parallel/ series implementatiuii . 

The following sections of the report concentrate: (1) on placing boundaries 

on the Instability in the paraliel/series implementation and (2) on evaluat- 
ing the stabili^ing effects of various system co. figuration changes, keeping 
in mind the critical differences between the implementations which apparently 
introduced the insLabil i ty . 


Effect of Input Amplitude and Bandwidth 

Addressed in this section of the report is the theoretical prediction 
of the stability boundary for sinusoidal and square wave Inputs, and the 
accuracy of that prediction, as verified by experiment. Using the iteration 
technique for describing function gain value as a function of system input 
discussed previously, the theoretical stability boundary as a function of 
input amplitude and frequency has been determined. This theoretical boundary 
is based on the value of bj /A for neutral stability determined from figure 
13# The experimental bou-'dary was found, as shown In the time histories, 
from the digital model. 

Sinusoidal Input Figure 23 illustrates the stability boundary for 
sinusoidal inputs of various magnitudes at the nominal limiter authority of 
20% (±.12 inches)# The input magnitude whicn the system will tolerate was 
found to be inversely proportional to the input frequency. From the figure 
it can be seen that the theory shows excellent correlation with experimental 
results for small input amplitudes, and becomes progressively less accurate 
as the amplitude increases. Recall from the previous comparison of predicted 
and actual responses that the theory consistently predicted a more stable 
response than experimental results showed, for the saturated limiter. 

Due to the distortion of the frequency response of ^ e nonlinear system 
between system input magnitude, M, and limiter Input magnitude. A, a dis- 
continuity or "jump” resonance, as described in reference 5, was found to 
exist. Shown in figure 24 is the frequency response of the linear system and 
the Jump discontinuity responses for three nonlinear system input magnitudes. 
At saturation (recall for nominal 20% series servo authority, the limits are 
±•12) for each system input magnitude, M, the limiter input magnitude. A, 
suddenly Jumps to a much higher value causing a sudden decrease in describing 
function gain, b^/A, and thus total system gain, sufficient to cause an 
unstable limit cycle. (Theoretically, this behavior has the characteristics 
of a hysteresis with the resonance frequency dependent upon the side from 
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which it is approached, but this was not successfully demonstrated experi- 
mentally.) Since the saturation point decreases with input magnitude, the 
resonance frequency also decreases as a function of input magnitude in a 
manner consistent with the inverse proportionality between tolerable input 
magnitude, M, and frequency, o), shown in figure 23. 

This type of insidious instability due to saturation would be extremely 
hazardous operationally because of its unpredictability to a pilot. Figures 
25 and 26 show the time history response which may be anticipated for inputs 
of 6i^|^ = ,4sin(1.22t) and .4sin(1.25t) , respectively. The difference 
between stability and instability in this case was found to be as little as 
.03 rad/sec. Because the transition to an unstable system is not gradual or 
predictable, in order to acquire an acceptable pilot-in-the-loop control 
system, the stability boundary mi. t be avoided during operation or eliminated 
entirely. 

Square Wave Input .— Sudden maneuvers require more, abrupt pilot inputs 
and are better approximated by a square wave. In this section of the report, 
the magnitude vs. frequency stability boundary established in the previous 
discussions has been extended to include a square wave. Using Fourier series, 
as c ascribed in reference 6, the square wave has been modeled analytically by 
equation (6) using up to the 21st harmonic of the fundamental frequency com- 
ponent . 

21 

6t * ~ ^ “ sin P t (6) 

'*’4; 7T " n 

n*l 

A theoretical model including the complexity of the 21st harmonic was 
chosen because of the very close resemblance of the time response using this 
representation and that of a square wave as shown in figures 27 and 28. In 
figure 28, upstream signals, acceleration, i|). Include a 

significant amount of the higher frequency Fourier series components. However 
yaw attitude response, is virtually identical to that commanded by the 
actual square wave, as in figure 27. 

System stability for a sine wave input was found to be very sensitive to 
input frequency. It would be anticipated that, because of the composition of 
very high frequencies required to model the straight sides of the square wave, 
system stability to a square wave input would be reduced significantly from 
that of a sine wave. 

Shown in figure 29 are the theoretical and experimental stability 
boundaries for a square wave system input. The theoretical curve was gene* 
rated using the Fourier series approximation as the input into the iteration 
technique for describing function gain value. The experimental curve reflects 
system response to an actual square wave, rather than the Fourier series 
approximation. As shown, the stability boundary is reduced significantly 
from the sinusoidal case in each magnitude ar.d frequency. For example, a 
square wave input of amplitude M/Mj^ « .1 becomes unstable at .15 rad/sec 
as compared with 1.8 rad/sec for a sinusoidal input of the same amplitude. 

As with the sinusoidal stability boundary, for higher input magnitudes, theory 
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predicts a more stable response than the experimental results show. Jump 
resonance phenomena also exist here, as shown in the distorted frequency 
response of system input magnitude to limiter input magnitude in figure 30. 
Due to the multi-frequency composition of the square wave, several jump dis- 
continuities exist at different resonance frequencies. 

Thus, the magnitude-frequency boundary is also a function of input 
shape; specifically, stability of the system response is significantly more 
sensitive to a square wave than to a sinusoidal input. 


Effect of Series Servo Authority 

One of the critical differences between the parallel/series and series 
servo systems is the higher limiter values of the latter which occur because 
of the proportionality between series servo authority and limiter magnitude. 
Nominal limits of the parallel/series system are based on 20% series servo 
authority and the maximum practical limits would be based on 50%. In this 
section of the report, the effect of variation of series servo limits on the 
nonlinear system has been evaluated. 

Theoretically, the result of increased series servo authority may be 
understood to have a stabilizing effect on system response by recalling the 
dependence of describing function gain, bj/A, on limiter size, a, from 
equation (3). As the limiter magnitude is increased, so is the minimum 
allowable describing function gain value, yielding greater system stability. 

From figure 31 (generated analytically) it may be seen that the sta- 
bility boundary as a function of series servo authority normalized to input 
magnitude varies almost linearly with low sinusoidal input frequencies. As 
the ratio of series servo authority to input magnitude increases (either by 
increasing series servo authority or decreasing input 'magnitude, as previously 
discussed) the system's tolerance to input frequency I increased. Note that 
a percentage of series servo authority normalized to input magnitude exists 
such that at or above which the system has a stable response at all input 
frequencies. 

TV » trend toward increased system stability for increased series servo 
authority is further Illustrated in figure 32. The comparison of magnitude- 
frequency tradeoff for nominal and Increased series servo authority, 20% and 
30%, respectively, was found to show a significant increase in the stability 
boundary at the maximum practical limit of 50% series servo authority. 


Effect of System Bandwidth 

Evaluated in this part of the report are the influences of system band- 
width on nonlinear system stability. Contributions to bandwidth that are 
considered are the yaw angle and rate feedback gains and the parallel servo 
gain in the forward loop. 
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Yaw Angle and Rate Feedback Gain Variations .— Variation of the yaw angle 
feedback gain, (as in figure 4), with simultaneous variation of the yaw 
rate feedback gain, , in order to maintain adequate closed loop damping, 
showed a significant increase in linear system phase margin bur almost no 
change in gain margin. This is illustrated in figure 33 for variations in 
and from the nominal values of 4.0 and 4.0 to .5 and 1.45, respec- 

tively. 

However, nonlinear stability boundary changes were found to be unfavor- 
able for reduction of yaw angle and rate f ecdboClv ^ct xU£i in figure 34. 

From the figure, two trends were noted: (1) as system bandwidth was reduced, 

the system's tolerance to magnitude at high input frequencies was increased 
slightly because of the lower describing function gain value for neutral 
stability due to the lower feedback gains, and K^. However (2), for the 
same amount of system bandwidth reduction, a much greater reduction in input 
bandwidth occurred. Equation (5) shows the relationship between system input, 
6iw^, describing function gain, bj/A, and the yaw angle and rate feedback 
gains, and K^. 

More favorable nonlinear system results were obtained by varying only 
the yaw angle feedback gain, while maintaining the yaw rate feedback gain 
at its nominal value. (This resulted in a 20% to 30% decrease in closed loop 
damping ratio.) Varying showed little effect on linear system stability 

criteria, as illustrated in figure 35 for values of from the nominal 

value of 4.0 to the limiting value of 0. 

Changes in the nonlinear system stability boundary with were found 

to exhibit the same two general trends as previously noted (figure 36), 
Specifically, these were a reduction in tolerable input bandwidth and an 
increased level of system tolerance to magnitude at high input frequencies. 

In this case, however, the increased tolerance to magnitude at high fre- 
quencies outweighs the input bandwidth loss, resulting in a net favorable 
effect on the nonlinear stability boundary. 

Recall that a critical difference between the parallel/series and series 
only systems was the introduction of an open loop pole at the origin into the 
parallel/series system by the parallel servo. This extra open loop integrator 
had a destabilizing influence on the parallel/series system as was shown in 
figures 13 and 14. Elimination of the yaw angle feedback (K^ » 0) changed 
the parallel/series system into one with only three open loop poles at the 
origin (as in the series only system). 

Figure 37 illustrates the variation of the system eigenvalues with 
describing function gain value for the closed loop system without yaw angle 
feedback. From the figure it may be seen that as the describing function 
gain, bi/A, approaches zero, the branches of closed loop eigenvalues depart 
vertically and remain in the left half plane. In figure 36, then, the sta- 
bility boundary is completely eliminated for * 0* Stabilizing the system 
in this manner does not preclude the existence of the low frequency nonlinear 
mode. The frequency and damping of this mode vary, as before, with describ- 
ing function gain value, calculated from equation (5) as a function of system 
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input. Eliminating the yaw angle feedba^'k and thereby stabilizing the system 
ensures only that the damping of this mode will be positive. 

Yaw angle feedback elimination changes the system from the designed 
rate-coiranand-attitude-hold configuration to rate-command only. System 
response, although stable for all inputs, lacks directional-hold capability. 
With an inadequate yaw rate sensor, heading drift may present a problem; also, 
removal of attitude-hold removes the control system^s resilience to wind or 
turbulence upsets. 


Parallel Servo Gaia Variations .— In terms of operational safety considera- 
tions, it was desirable to maximize the saturation time following a runaway 
failure, and thus minimize parallel servo gain, Kp . Another consideration 
in gain selection was the upstream limiting of the parallel servo input signal 
by the series servo lijniLer. Due to the proportionality of the series servo 
limiter value to series serve authority, the parallel servo gain was chosen 
to be inversely proportional to this parameter. Parallel servo gain, then, 
ai: shown in equation (7) is a function of the time which is allowed for satu- 
ration of the parallel servo in the case of a runaway failure and of the 
fraction of the total control power allotted to the series servo. 


^parallel 

servo 


1 


/ Series Servo\ 

/Saturation\ 

I Authority ) 

1 Time ) 


(7) 


The same trends as in the previous discussion were apparent with the 
system bandwidth reduction due to lowering of the parallel servo gain, Kpg. 
Linear system phase and gain margins were found to increase, although much 
more significantly than in the case of the attitude feedback gain, with 
parallel servo gain reduction, as illustrated in figure 38. From the figure, 
it may be seen that a significant stabilization of the linear system occurs 
with reduced parallel servo gain. 


With the nominal value of 20Z series servo authority, the minimum value 
of parallel servo gain may be established to be Kpg * .25, based on a maxi- 
mum practical saturation time of 20 seconds. If the series servo authority 
were increased to the maximum reasonable level of 50Z, a minimum value of 

is acceptable. The same general trends, specifically a reduction 
in tolerable input bandwidth and an increased level of system tolerance to 
magnitude at high input frequencies as were previously noted for attitude 
feedbo. gain reduction, were found to also apply to the reduction of parallel 
sf^/o gain. These trends are Illustrated for the nonlinear system with 20Z 
,erles servo authority in figure 39 where it may be seen that the undesirable 
loss of input bandwidth far outweighs the gain in system tolerance to magni- 
tude at high frequencies for a given reduction in parallel servo gain. 

Figure 40 shows the change in nonlinear etability boundary with parallel 
servo gain for a system with 50Z series servo authority. The combination of 
the Inherently stabilizing effect of Increased servo authority and the smaller 
minlmt'm acceptable parallel servo gain of Kp^ * .1 resulted in a net favor- 
able effect on the nonlinear stability boundary. 
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CONCLUSIONS 


Tht. analysis described in this report encompassed linear and nonlinear 
stability analyses of the state-rate-feedback-implicit^model-following control 
system of reference 3. Theoretical stability predictions of the nonlinear 
system were accomplished using a describing function linearization technique. 
F"<^perimental results obtained from a digital model were used to substantiate 
the theory. From this analysis, the following specific conclusions have been 
drawn: 


1. The experimental analysis using the simplified Fortran model was 
accomplished non-real-time with a sampling period of T * .005 seconds. 
Although the real-time computer simulation uses a sampling period of T = .05 
seconds, the linear system’s dynamic characteristics were found to be repre- 
sented much more accurately with T = .005 seconds. 

2. Analytical prediction of nonlinear system stability using the 
describing function linearization technique for a simple limiter was found to 
be valid for each of the sinusoidal and square wave system inputs. 

3. In its nominal configuration, the full authority series servo mode 
was found to be stable at all system inputs; the parallel/series servo mode 
was found to be conditionally stable and possibly unacceptable for inputs 
which saturate the limiter upstream of the parallel and series servos. This 
conditional stability results from the extra open loop pole at the origin 
introduced in the parallel/series implementation. The effect of the series 
servo position limiter is sufficient to reduce the system gain and introduce 
a low frequency Instability. 

4. A stability boundary, in which tolerable input magnitude was 
inversely proportional to input frequency was found to exist in the 
parallel/series servo mode. In terms of this stability boundary, signifi- 
cantly more sensitivity was observed to a square wave input due to its high 
frequency composition than to a sine wave input. 

5. Increasing the fraction of total control power allotted to the 
series servo was found to have a significantly favorable effect on the sta- 
bility of the nonlinear system with parallel/series servo implementation. 
Increasing the series servo authority to its maximum practical level was 
found to raise, but not eliminate, the stability boundary. 

6. In terms of the magnitude-frequency stability boundary, reduction of 
yaw angle and rate feedback gains was found to affect the nonlinear stability 
boundary unfavorably. Yaw angle feedback gain reduction (alone) Improved 
stability, resulting in a relatively small input bandwidth reduction and a 
more significant increase in system tolerance to magnitude at high input fre- 
quencies* Elimination of yaw angle feedback was found to stabilize the 
parallel/series system; however removal of the attitude hold capability 
reduces the system's resilience to rate gyro drift or atmospheric distur- 
bances . 
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7. Parallel servo gain reduction (alone) was found to have an adverse 
effect on nonlinear system stability. The small increase in system tolerance 
to magnitude at high input frequencies was by far outweighed by the rela- 
tively large loss of input bandwidth. The combination of parallel servo gain 
reduction and series servo authority increase was found to improve nonlinear 
system stability. 
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APPENDIX 


GAIN EQUIVALENCES AND INITIAL VALUES 


SRFIMF Initial Numerical Value 


KlO 

K20 

*^30 

FCN 


± 2.1 in. 
0.3948 
0.9349 
4.0 

4.0 
0.3324 

± 0.12 in. 
0.05/sec 
2.533 

1.0 
1.0 
7.5 


Analysis 

SRFIMF 

Initial 

Numerical Value 



— 



±2.1 in. 

Kl 

0^ ^ k^h> 

1.33 

K 2 

(Kio + K 20 ) 

1.33 

K 3 

(K 30 X X FCN X KipoRQUE^^z^ 

0.7082 


^*^ORQUE^^z^ 

0.0373 

*S>8 

— 

2.5 


(K 30 X Kpg X Kyc X FCN) 

47.49 

*^88 


7.0 

Ks 

(Kgs K 30 ) 

17.73 



4.0 


H 

4.0 
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Analysis 

SRFIMF 

Initial 

Numerical Value 

±a 

±Uitl>c 

±0.12 In. 

^7 

"7 

0.05 /sec 

^2 


0.10/ sec 

^torque 

— 

1200 deg-lbf 



32192 slug-ft^ 



TABLE 1.- C(»(PARISON OF THEORETICALLY PREDICTEL AND EXPERIMENTAL SYSTEM RESPONSE; 
NOMINAL 20% SERIES SERVO AUTHORITY, a - ±0.12. 


































PEDALS 



Figure 3.- Simplified yaw control system, mechanization and actuation dynamics 

























Figure 6.- Root locus of linear i^ystem^ parallel/series servo mode 
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Figure 7,- Partial root locus of linear system, parallel/ seri zs ser 'o mode. 
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Figure 11.- Time history vs. linear analysis, parallel/series servo mode. 
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Figure 13.- Variation of closed- loop roots with magnitude of nonlinearity (b./A), 

perallel/serles servo mode. 
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Figure 16.- Nominal configuration, parallel/sariaa aarvo, 6j • 0.4 aln (l.Ot). 
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Figure 19.- Nominal configuration, series servo, 5. ~ 0.1 sin (l.Ot). 
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Figure 20.- Nominal configuration, series servo, = O.A sin (l.Ot). 
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Figure 21.- Nominal configuration, series servo, 6 =0.1 sin (8.0t). 
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Figure 23.- Amplifide-f requency stability boundary for .sinusoidal Input. 
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M = 0.4 
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(jj, rad/sec 
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(INCLUDING HYSTERESIS) 


M = 0.2 
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O CL 
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SYSTEM 
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i^igULe 24.- Frequency response of linear and nonlinear systems; inset shows Jump 
discontinuity upon saturation; system input ~ * M sin ut; 

20% serjps servo .lui!:.-!. Ity; ±0.12 in. limiter '^authority. 
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Figure 25.- Nominal configuration, parallel/series servo, 6_ ■ 0.4 sin (1.22t). 
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%AX = 2.1in. 

20% SERIES SERVO AUTHORITY 


±0.12 in. LIMITER AUTHORITY 
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Figure 29.- Amplitude- frequency stability boundary for square wave input. 
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Figure 30.- Frequency response cf linear and nonlinear systems, including jump 
resonance discontlruity upon saturation^ 
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Figure 32.- Amplitude-frequency stability boundary as a function of series servo 

authority. 
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Figure 
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Figure 35.- 


car system phase ^ j 
feedback gain, K , " ®3rglns with yai 



5i =M sin cot 
%AX " in. 

20% SERIES SERVO AUTHORITY 
±0.12 in. LIMITER AUTHORITY 



u), rad/sec 

Note: With =» 0 => stability boundary eliminated. 


Figure 36.- Amplitude-frequency stability boundary as a function of yaw angle 

feedback gain, K^. 
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Figure 38.- Variation of linear system phase and gain margins with parallel servo 

galnt I^ps* 
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